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Abstract 

In this pedagogical work we point out a subtle mistake that can be done 
by undergraduate or graduate students in the computation of the electrostatic 
energy of a system containing charges and perfect conductors if they naively 
use the image method. Specifically, we show that naive expressions for the 
electrostatic energy for these systems obtained directly from the image method 
are wrong by a factor 1/2. We start our discussion with well known examples, 
namely, point charge-perfectly conducting wall and point charge-perfectly con- 
ducting sphere and then proceed to the demonstration of general results, valid 
for conductors of arbitrary shapes. 
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1 Introduction 



The typical problem in electrostatics consists of determining in all relevant space the 
electrostatic field generated by some set of charges from the known charge distribution 
itself as well as from the appropriate boundary conditions pertaining to the situation. 
Solving this problem amounts to finding a static potential $(x) obeying Poisson's 
equation, namely, 



subject to the suitable boundary conditions. Although one can approach this task in 
many different ways - whose convenience depends on the particular problem being 
dealt with -, Poisson's equation (given definite boundary conditions) has a unique 
solution for each charge distribution p(x). This allows one to look for solutions in 
any desired fashion: if one finds a potential that obeys both the boundary conditions 
and Poisson's equation for the correct p(x), it must be the correct potential for the 
given configuration. 

A especially suitable method for simple situations with point charges and dipoles in 
the presence of conductors is called the image method. It consists of finding a different 
configuration in which the conductors are replaced by some charge distribution so 
that the potential of the entire setup in the region of physical interest (i.e., outside 
the conductors) fits the appropriate boundary conditions and also obeys Poisson's 
equation. These fictitious charges put in place of the conductors are called image 
charges. Since any potential created by charges obeys Poisson's equation (eqil]), once 
the boundary conditions are satisfied, the field created by all charges (real and image) 
is the field obtained in the actual configuration outside the conductors. It is important 
to state that in the fictitious configuration the charges outside the conductors must 
be in the same place as in the real distribution, or else one would find a solution 
to Poisson's equation with a different p(x): all image charges must be placed in the 
space originally occupied by the conductors. 

The method can then provide the force acting on a charge in the presence of 
conductors. However, if one tries to naively use the image method to compute the 
electrostatic energy of that configuration, an incorrect result will be found, namely, 
one arrives at twice the correct energy, as Griffiths [Ij and Franklin [2J have shown for 
particular cases. In simple situations as the point charge-plane wall setup, it is rather 
easy to realize that the real configuration has half the energy of the fictitious one, but 
in more involved and less symmetrical geometries this is no longer an obvious issue. 
In the case of a charge in front of a sphere one could lucidly expect to find a prefactor 
depending on the sphere radius which would only tend to 1/2 if the radius tended to 
infinity (reobtaining the charge- wall result). 

Our goal is to show that, whatever the shape and amount of perfectly conducting 
bodies near a point charge, the electrostatic energy of the system is 1/2 times the 
Coulombic energy of the interaction between the real charge and each image charge 
of the problem. We shall then generalize that result for more than one source charge. 

Our article is organized as follows: we begin in Section 2 with the familiar problem 
of a charge and a conducting wall, and also comment on the case of a wedge. We then 
proceed, in Section 3, to a less symmetrical geometry, that of a conducting sphere. 
Section 4 is dedicated to the theorem demonstration in the general case, and we leave 





the last section for conclusion and remarks. 



2 Usual case: point charge and conducting wall 

In order to state the problem and emphasize the important point in the clearest way 
we start our discussion by considering in this section the simplest problem that can be 
solved by the image method, namely, to find the electrostatic field of a point charge 
in the presence of a perfectly conducting wall. This problem can be found in many 
standard textbooks [H El El H] , so we go directly to the point here. For convenience, 
let us choose our cartesian axes OXyZ such that the region z < is filled with a 
perfectly conducting material and a point charge q is at position (0,0,z), as shown 
in Fig. [D 
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Figure 1: Real and image charges for the charge- wall case 

According to the image method, the force exerted on the real charge q by the super- 
ficial charge distribution induced on the surface of the conductor is given by 

One could naively think that the electrostatic energy of the system shown in Fig. [1] 
would, too, be given simply by the Coulombic energy between point charge q and 
its image, — g, namely, U = — (g^/47reo)(l/22;). However, this is not true, as we can 
easily verify if we take the gradient of U, 

.^U^-.1(-J^^)^-2(JI-J^),^,,. (3) 
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The correct expression for U has an additional factor of 1/2 and can be readily 
obtained if we start by the very definition of U as the total external work to bring all 
the real charges (including the surface distribution) from infinity to the final static 
configuration Cf. The easiest way to compute it is to picture the setup of all charges of 
the system along their paths from infinity to C/ always consistent with the presence 
of the conductor. This guarantees that the surface charges move across regions of 
constant potential, without any work required to bring them. Under these conditions 
(in the electrostatic context radiation effects are negligible), we have 

U = -W,,t = - Fg{r')-dr'. (4) 



Substituting eq.(j2]) into the above equation, we get 



In other words, the electrostatic energy of the system formed by point charge q and 
the conducting region is half the electrostatic energy of a point charge q located at 
(0, 0, z) and a point charge —q located at (0, 0, —z). This result is in agreement with 
([2]) as can be readily seen. This kind of discussion can be found in many textbooks, 
like Griffiths's [I], among others. 

In this simple case, the factor 1/2 could also be anticipated by symmetry argu- 
ments, as follows. First, recall that (l/2)eoE^ is the energy density of the electrostatic 
field. With this in mind, we easily see that the energy of a system formed by the 
charges q at (0, 0, z) and —q at (0, 0, —z^ (with no conductor at all) is equally divided 
between the regions z > and 2; < 0. In this calculation, we must, of course, exclude 
the self-energies of each charge. 

Symmetry also allows us to deal with the case of a charge near an infinite wedge 
whose aperture angle equals vr/n for any positive integer n in = \ corresponding to 
the plane wall). In these cases, the entire space can be divided into In sectors with 
that same angle, one corresponding to the outside of the conductor, the remaining 
2n — 1, to the space filled by the conductor. An image charge will be in each sector, 
except for outside the conductor, where the real charge q\ lies. We shall label the 
sectors and the pertaining charges with integers, i = \ refers to the real charge, 
i = (2, ...,2n), to the image ones. Symmetry allows us to say that the configuration 
energy would be l/2n times the Coulombic energy of all 2n charges, i.e.. 



2n 2n ^ 2n 2n 
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where Vj (ri) is the potential created by charge qj at the position ri of charge and 
is the distance between charges qi and qj. This seems to indicate that the prefactor 
depends on n, i.e., on the wedge angle at hand, but we must pay closer attention 
to the expression we are comparing the energy with. The interaction energy of the 
charge- wedge system is l/2n times the energy of the total system composed of 2n 
charges, while we wish to compare it with the Coulombic interaction between the real 
charge and each of the images (image-image interactions not being included). This 
distinction can be very subtle; in the case n = 2 we wish to compare the actual energy 
with the energy of the pairs (51,52), {QijQs), {qi,Q'i), while the double summation on 
eq.([n]) includes these plus the pairs (52, 53), {q2, Qi), {is, Qi)- More generally, the double 
summation on eq.d^]) comprises n{2n—l) pairs. If one only counts the interaction 
between the real charge and each image, one finds (2n— 1) pairs. Moreover, using the 
fact that odd-numbered charges have the value q and even-numbered ones, the value 
—q, together with the symmetry of the configuration, one can see that 

2n 2n 2n 
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An interested reader may verify eq. ([7]) for any particular value of n. We thus conclude 
that 

^ 2n ^ 2n 

i=2 i=2 

and once more the energy of the configuration is half of the interaction energy between 
the real charge and each image. 

A natural question then arises: what happens in less symmetric situations or even 
in situations where there is no symmetry at all? From now on, this answer is our main 
concern. However, we shall do that in two steps. First, we shall consider in the next 
section another example and work out the result explicitly. Then, we shall attack a 
completely general situation of one (or more) charged particle (s) in the vicinity of 
grounded or neutral perfect conductors of arbitrary shapes. 



3 Point charge and a perfectly conducting sphere 

Let us consider as our next example a point charge q near a perfectly conducting 
grounded sphere of radius R. Suppose the distance from charge q to the center of the 
sphere is a, a > R. For simplicity, we choose the axis OX with its origin at the center 
of the sphere so that the position of charge q is given by (a, 0, 0), as shown in Fig. O 
It is well known that the surface charge distribution on the sphere is such that the 
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Figure 2: Real and image charges for the charge-sphere case. The image charges are 
q' = —qR/a and —q'. 

force on q is the same as if there were no sphere at all and a charge g' = —{R/a)q 
were located at {R^/a, 0, 0) (see Ref. [I]). Therefore, the force exerted by the surface 
distribution of the sphere on the point charge q is given by 

^ q{-qR/a) x ^ -q^R ax 
' 47reo (a-i?Va)2 ineo {a? - R'f' ^ ' 



According to the previously presented discussion, we do not expect the electro- 
static energy of the charge-sphere system shown in Fig. [2]to be given by the Coulombic 
interaction energy between real and image charges. However, we have no reason, a 
priori, to say that the correct answer is obtained simply by including an additional 
factor 1/2 as occurred in the cases discussed in Section 2. It would be natural, though, 
to expect a factor depending on R and a that, for R —>■ oo with R — a kept constant, 
reduces to the previous factor 1/2, since that limit reproduces the charge- wall case. 

Let us then perform the explicit calculation using, as before, the very definition 
of the electrostatic energy of a configuration in light of the comments made before 
eq.(jl]). Doing that along the axis OX and using eq.(l9]), we may write 

u = -Wext = r dx 

47reo L {x' - R^f 

1 q{-Rq/a) 

2 47reo (a - i?Va) ' ^ ^ 

which is nothing but 1/2 the Coulombic energy between charge q and its image 
q' = —Rq/a. At first sight, it seems amazing that the same factor 1/2 appears. 
This suggests that this will happen for conductors of general shapes. In fact, this is 
precisely what happens, as we shall demonstrate in the next section. 



4 General case 
4.1 One source charge 

We shall now consider one point charge q in the vicinity of a set of N perfect con- 
ductors of arbitrary shapes, which can be either neutral or grounded. Let xq be its 
position in space with respect to some reference frame. Our purpose here is to ob- 
tain an expression for the electrostatic energy of this configuration in terms of the 
Coulombic interaction energy between q and each image charge necessary to solve the 
problem, which would be the energy necessary to bring in the charge q from infinity 
with every image held fixed at its final position. For convenience, we shall anticipate 
the final result in the form of a theorem, namely. 

The electrostatic energy of a point charge q near N perfect conductors of 
arbitrary shapes, each conductor being either neutral or grounded, is half 
the Coulombic energy between the charge q and each image charge. 

We now present a simple demonstration of this theorem. It is convenient to 
start by the following expression for the electrostatic energy [/ of a general charge 
distribution, 

f/=^^p(x)$(x)rf\/, (11) 

where p(x) is the charge volumar density at position x, $(x) is the electrostatic poten- 
tial at position x due to all charge distribution and 7^ is a region of space containing 
all charges. Of course, whenever point charges are present in the distribution, we 
must subtract from the above expression the corresponding infinite self-energies. It 



has already been shown in the hterature (for instance, in [T]) that this expression is 
equivalent to computing the external work to bring the charges from infinity. 

Since our distribution consists of a point charge q located at position xq and 
surface charge distributions on the conductors, eq. (|TTi) can be written as 

1 I ^ f 

-q^^,) + -J2f afc(x)$(x)t/Afc, (12) 
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where $(xo) is the potential at position xq created by all charges of the system 
except point charge g, describes the charge distribution on the surface of the 
k-th conductor. Writing $(xo) instead of $(xo) in the first term of the r.h.s. of the 
previous equation is equivalent to subtracting the (infinite) self-energy of the point 
charge q. Recalling that each surface Sk is an equipotential surface, whose potential 
we denote by we obtain 



1 1 ^ 

U=-q <l>(xo) + ^Y^^kf a,.(x) dAk. (13) 

k=i ^'^k 



Since each conductor is either grounded, from which $fc = 0, or neutral, from which 
the surface integral is zero, the last term on the r.h.s. of eq. (fT3|l always vanishes, and 

f/=lg<l(xo) (14) 

Now, all we need to do is to invoke the image method to finish our demonstration. 
Image charges are, by definition, imaginary charges situated in the nonphysical regions 
(inside the conductors) that create at any point of the physical region (outside the 
conductors) the same field as created by all surface distributions of all conductors. 
Hence, $(xo) is precisely the electrostatic potential at position xq due to all image 
charges, so that we can write symbolically, 

U = ^q <i>(xo) = \u(^q; {images}) , (15) 

where u(^q; {images}) means all Coulombic interactions between point charge q and 
each image. This completes the demonstration for one point charge q. 



4.2 Many source charges 

The result can be further generalized to accommodate the presence of more source 
charges. We shall now consider a set of M point charges qi, q2, qM in the vicinity 
of a set of N perfect conductors of arbitrary shapes, which can be either neutral or 
grounded. Let Xj be the position of charge qi in space with respect to some reference 
frame. The expression for the electrostatic energy of this configuration shall also 
include terms due to interaction of the real charges. Anticipating the final result once 
more: 

The electrostatic energy of a set of M point charges qi,q2, ■■■,qM near N 
perfect conductors of arbitrary shapes, each conductor being either neutral 
or grounded, is the Coulombic interaction energy between the real point 
charges plus half the sum, from i = l to i = M, of the Coulombic energies 
between charge q^ and each image charge. 



Let us demonstrate this second theorem. We start once again by eq. lfTTi) . 



U=l f p(x)$(x)rf\/. (16) 

We can write the expression separating p(x) in the contributions due to each point 
charge located at position Xj and due to the surface charge distributions on the 
conductors, 

M N 

i=i k=i 

While (Tfc describes the charge distribution on the surface of the k-th conductor, 
we subtract the self-energy of point charge gj by substituting the potential $(xj) by 
^^(xj), the potential at position Xj created by all charges of the system except point 
charge itself. The second term in the r.h.s. of this equation vanishes analogously 
to how the second term on the r.h.s. of eq. (fT2!) does. We obtain 

1 *^ 

u = -Y,q:M^i) (18) 

2=1 

The potential can be split in the potential due to the other point charges qj and 
the potential due to all surface charges together, 

MM M 

f'-^EEi^^^ + iE^Wx.). (19) 

i=i j=i * J I i=i 

The first term can be readily recognized as the Coulombic energy between the source 
charges, as can be seen in the literature fl] (we remember that each pair is counted 
twice in that double summation). Invoking the image method as before, we can state 
that $5u^.j(xj) is the electrostatic potential at position Xj due to all image charges, 
allowing us to write, symbolically, 

1 

U = ?7^{all source charges} j H — u(^qi; {images}^ , (20) 

i=l 

thus completing our most general demonstration. 

Eq. ( l20l) can be interpreted in terms of pairwise Coulombic energies as follows: the 
energy of each pair composed of two real charges enters the expression integrally; 
pairs that comprise a real charge and an image one take on a factor one half; pairs of 
image charges are not present in eq. (l20|) . 



5 Conclusions and final remarks 

We have calculated the electrostatic energy of systems composed of a point charge 
and conductors of various geometries. We started with the simple case of a plane wall 
and obtained as a result one half of the Coulombic energy between real and image 
charges, which could be easily understood considering the fields' energy density and 



the symmetry of the problem. We then indicated, solely on symmetry arguments, 
that the energy of a system composed of a charge and a wedge of aperture angle 
7r/n [n a positive integer) should also be one half of the Coulombic energy between 
the charge and every image. We then moved on to a less symmetrical geometry, 
a spherical one, and there, too, we found the same factor 1/2 when comparing the 
energy of the actual system to the pairwise Coulombic energy between the real charge 
and each image. We then proved the theorem that this same factor one half arises 
in every problem of a point charge in the presence of conductors that is solvable by 
the image method, whatever its geometry may be. The argument was completed by 
generalizing the theorem for the case of more than one source charge. 

In many situations involving more than one source charge, the interest lies on the 
interaction energy between the sources and the conductors, but not on the self-energy 
of the sources. For those cases eq. (!20|) is best written in the form 



which suits multipole sources especially well. 

Although our discussion was in the classical context, it may serve as a use- 
ful guideline for quantum problems involving dispersive forces between polarizable 
atoms/molecules and conducting bodies. 
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